Lattice four-dimensional A/=4 SYM is practical 
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We show that nonpertubative lattice studies of four-dimensional Af—4 Super- Yang-Mills are within 
reach. We use Ginsparg- Wilson fermions to avoid gluino masses and an exact implementation of the 
(chiral) /^-symmetry, which greatly limits the number of counterterms that must be fine-tuned. Only 
bosonic operators require fine tuning, so all tunings can be done "offline" by a Ferrenberg-Swendsen 
type reweighting. We show what measurables can be used to perform the tuning. 

PACS numbers: 



It is often assumed that the study of supersymmetry 
(SUSY) on the lattice will be difficult in four dimensions 
because multiple relevant/marginal operators that vio- 
late SUSY are allowed by lattice symmetries. A notable 
exception is pure Af=l super- Yang-Mills (SYM). In that 
case, lattice chiral symmetry in the form of Ginsparg- 
Wilson (GW) fermions prevents additive renormalization 
of the gluino mass in the continuum limit, and hence the 
only SUSY-violating operator is forbidden in that limit 
by setting the bare mass to zero [3 N=A SYM has scalars 
and Yukawa couplings, and so various masses and cou- 
plings will receive divergent corrections in the continuum 
limit and must be (nonperturbatively) fine-tuned away. 
Hence lattice Af—4 SYM has always seemed impractical 
by this fine-tuning approach. [l3j 

We show here that this is not the case. Using GW 
fermions the four gluinos can be kept massless and the 
SU(A) R symmetry in the lattice theory can be preserved, 
which greatly restricts the form of the renormalizations. 
This leaves a scalar mass, two quartic couplings, and a 
Yukawa coupling to tune. The Yukawa coupling can be 
tuned by rescaling the scalar kinetic term (see below); 
so all tunings can be done by adjusting bosonic terms in 
the action. This allows the tunings to be done by the 
"Ferrenberg-Swendsen method" 0, exploring a wide 
swath of coupling constant space "offline" from the re- 
sults of a single Monte-Carlo simulation. The parameter 
range available with good statistics can be enlarged using 
multicanonical techniques [4j, [5| . Thus we arrive at the 
encouraging result that all fine-tuning can be performed 
through an "offline" analysis, i.e., new simulations and 
fermion matrix inversions are not required. 

The continuum field content is SU(N C ) YM theory 
with an SU(A) R internal symmetry; there are four Majo- 
rana fermions whose left handed components transform 
in the fundamental 4 representation of SU(4) R and 6 
real scalars in the antisymmetric tensor representation. 
The six real scalars will be expressed with a single index 
<j) m , m—1. . .6, or composed into SU(4:) R Weyl matrices: 
&i 3 =(t>rnV m ,i 3 and y =(/) m c7^, where ct's are just SU(4) R 
Clebsch-Gordon coefficients involved in 4*3 6®4. The 



continuum action is 
S = ±Tr{±F 2 + \DM 



(1) 



+V2^ (<j> ij PL-(<p ij yP R ) i>j + [^m,0«Pm,0n]}- 

The SU(4)~SO(6) preserving bosonic lattice action is 
a trivial transcription from the continuum. Of course we 
must allow for generic coefficients and non-SUSY terms, 
so that the SUSY-restoring counterterms can be tuned. 
In our case these are entirely scalar terms. 

It has been argued that naive lattice Yukawa terms 
lead to inconsistencies in either the chiral or Majorana 
projections (dependingon how the Yukawas are tran- 
scribed to the lattice) [6|], so the fermionic implementa- 
tion is more subtle. Following Liischer [7] and Kikukawa 
and Suzuki Q, we introduce a set of auxiliary fermionic 
fields 'F, with lattice fermionic action 

S F = a 4 Tr{V^ Gw ^ - a" 1 * (2) 
+yV2(^+f )t (<j>«P L - (^yP R ) ftM-*)i}. 

This action possesses an exact SU(A) R symmetry, with 
the scalars transforming as in the continuum and the 
fermions transforming according to 



(ty/ie = (TP L -T*P R )^ , 

SV/ie = (T+T*) 75 D GW ^ + (T Pl—T* Pr)^ , 

= i> (T*P L -TP R ) + *(T+T*) 75 , 
5^ fie = -<F (TP L -T*P R ) . 



(3) 



Here P, 



L/R 



±(1±7 5 ) = i(l±7 5 (l-2D QW )) are the lat- 
tice modified chiral projection operators, T is the gener- 
ator of SU (4)r in the fundamental (4) and we have sup- 
pressed the SU(A) R indices. Hence (V'+'F) and (yjj+^f) 
transform like the continuum ip, vb fields. 

This auxiliary field method preserves the i?-symmetry 
exactly and keeps the Yukawa terms ultralocal. It is 
also consistent with the Majorana decomposition, so 
the fermionic determinant is an exact square; taking its 
square root to implement the Majorana nature of the 
fermions retains locality. The cost is the introduction of 
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an extra fermionic excitation "J, which is however nondy- 
namical with C(a _1 ) mass, so it decouples from the the- 
ory in the continuum limit. Note however that in building 
lattice operators which contain fcrmions, we must always 
use the combination tp+^l, since this is the combination 
with correct i?-symmetry transformations. 

We are interested in tuning the lattice action such that 
the effective infrared description is jV=4 SYM. Generi- 
cally there is a nontrivial matching between lattice and 
effective IR theories and all relevant or marginal terms 
consistent with lattice symmetries will appear in the in- 
frared, except at special points in bare parameter space. 
We can arrive at the desired special point (i.e., 7V=4 
SYM) by introducing the SUSY-violating operators into 
the bare action and fine-tuning counterterms. These 
counterterms fall into three categories: a scalar mass 
term, a Yukawa term, and two or four scalar quartic 
terms, depending on the number of colors for the gauge 
group, restricted here to SU (N c ). In particular, if N c < 3 
then SU(&)r symmetry allows only two unique quartic 
terms: (Trcj) m (j) m ) 2 and (Tr(/) m ^„) 2 (trace over gauge in- 
dices, R- indices explicit). For N c > 3 the quartic invari- 
ants Tv(j) m (j) m ip n (l) n and Tr<fi m (j) n (f) m (j) n are also indepen- 
dent (for N c < 3 they are linearly dependent on the first 
two). 

Rescaling the Yukawa term, schematically yipcjnp, can 
be accomplished through a rescaling of the scalar ki- 
netic term. Consider replacing yipifnp + \Dcj)\ 2 with 
yipipip + Z\D(j)\ 2 . In terms of the canonically normal- 
ized scalar field this is (y / ' y/Z)jp(j)ip + \D(f>\ 2 . Therefore 
we include a counterterm for the kinetic term rather than 
for the Yukawa. (Note that the rescaling of the scalar ki- 
netic term also rescales the scalar potential; but there are 
distinct counterterms to undo this rescaling separately). 

Suppose we perform a Monte Carlo simulation at one 
value mi of the scalar mass m, so that the configu- 
rations sample the distribution determined by the ac- 
tion S = S m= o+5 Jm 2 <j) 2 . Following the "Ferrenberg- 
Swendsen reweighting" method 0, (3[ one can use the 
following "reweighting identity" to compute the expecta- 
tion value of an operator O for the distribution with a 
mass 7712 : 

There is a limited regime of utility to this technique, 
due to the so-called "overlap problem." For instance, if 
the exponential in ([J} is large where the simulated distri- 
bution has little weight, a finite sampling will have large 
errors. The mismatch of the distributions gets worse as 
the number of lattice sites increases, because the expo- 
nent is extensive. 

A way to ameliorate the overlap problem, which has 
been found to work in other contexts, is "multicanonical 
reweighting" One replaces 5" with S + W[Oi, 2 , • • ■], 
where W[Oi, O2, ■ ■ •] is a carefully chosen function of 



some small set of observables (in our case W will be a 
function of / (j) 2 , the distinct / <^ 4 's, and J(D(j)) 2 ). The ex- 
pectation value of an observable in the distribution cor- 
responding to S is: 

{0)=J2 c Oce w ^/j2 c e Wl0 ^- ] (5) 

W[..] produces a weighted average over a continuum of 
canonical ensembles, some of which will have a good over- 
lap with the distribution that one is reweighting to. The 
challenge is to design a W such that sampling is flattened 
over the range of observables one is interested in. 

Two approaches to engineering a good function W ex- 
ist: (1) a bootstrap method that iterates between Monte 
Carlo simulation and adjusting W, and (2) optimizing 
W w.r.t. its parameters, in a small volume, and then 
using step-scaling to extrapolate to a good estimate for 
W in larger volumes. For instance one can start with 
4 4 and 6 4 volumes, where statistics accumulate rapidly 
and unreweighted simulations still cover broad parameter 
ranges. 

In our case the reweighting function W will depend 
on the four bosonic contributions to the action, /<?!> 2 , 
f(D(f)) 2 , J and J (the two quadratic and two inde- 
pendent quartic operators you can form from the scalars, 
integrated over space). By sampling with the weight 

p = T>[A,<j>]det[D] e S[A,*;»%,-] e -w[f4> 2 ,...] j ^ 

and you can reproduce the ensemble at some particular 
set of values for m 2 , Z*, A and A', via 

Z = J2 c e +w[ -- ] e -[( m2 - m o)^+(^-^)(D0)^+...] ^ ^ 

with W chosen so that the sample has a reasonable num- 
ber of configurations for all values of Jcf> 2 , J(D<j)) 2 , J cf) 4 
and /(</> 2 ) 2 within some interesting range. 

Now, define the gauge invariant effective potential in 
finite volume as follows: 

e - QV [A*] = j2 c e+ w ^ e-U m2 - m ^-^xS(A 2 -^) (8) 

with il the 4-volume (so that A 2 represents the mean 
value of the squared scalar field) . It is easy to check that 
varying m 2 changes V [A 2 ] by adding a linear component. 
Therefore, measuring T^[A 2 ] immediately determines how 
(<j) 2 ) and F=lnZ vary as a function of m 2 . One easily 
generalizes to the effective potential as a function of all 
four bosonic operators, which gives a quick way to explore 
the effect of varying parameters. 

It remains to specify how to tune the parameters. Con- 
sider A/"=4 SYM with an added a 2 (Tr (/) m (j> m ) 3 term in 
the potential, but deformed by mass and quartic interac- 
tions. The SUSY point is a second order phase transition 
point; for negative quartic deformation there is a first or- 
der transition as m 2 is varied (the system jumps from a 
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massive state about the origin to a massive state with a 
large VEV because the quartic term bends down). For 
positive quartic deformation there is no phase transition. 
The optimal value of m 2 (fixing other parameters) is the 
point of maximal susceptibility ((J^ </> 2 ) 2 ) — ((Xc^ 2 )) 2 - 
Determining this point in finite volume leads to 0(1/ L 2 ) 
errors in the determined value of to 2 , which can be im- 
proved by scaling over multiple volumes. Finding the flat 
quartic term which gives second-order behavior should 
also be possible; it has been successfully achieved in the 
context of the electroweak phase transition . Therefore 
it should be possible to use the phase diagram to tune 
at least two parameters. Note that we needed to add a 
lattice-size 4> 6 term; this is harmless but it raises the issue 
that the flat direction actually means that unbroken 7V=4 
SYM is not well behaved in finite volume; the moduli are 
not fixed and the partition function diverges because of 
the integral over the infinite moduli space. Therefore it 
will always be necessary to break SUSY somehow. We 
advocate doing so via twisted boundary conditions; for 
instance, instead of periodic boundary conditions we can 
add a rotation by angle O to all fermionic fields in one 
direction. The choice 0=7r is the maximal global break- 
ing of SUSY and corresponds to treating the thermal en- 
semble; intermediate values of O break SUSY by smaller 
amounts. This lifts the moduli degeneracy without any 
local SUSY breaking; the effects of O are only visible in 
correlations at the scale of the lattice size, which is any- 
way contaminated by being in finite volume. 

For the case N c > 3 the effective potential should show 
multiple flat directions in the space of quartic operators; 
only one quartic direction (some linear combination of 
the input quartics, due to mixing) should rise steeply. 
Therefore we expect it should be possible to tune the 
"extra" quartic operators in the case N c > 3, leaving 
only one quartic and the Yukawa coupling/wave function 
to tune. 

If SUSY is exact then the (i?-symmetry 4) supercur- 
rent S^j is conserved, so (d^ S^i(x)0(y)) vanishes at 
Xy^y for all local operators O. We can use this property 
to measure whether we are at the SUSY point in param- 
eter space, and therefore to tune parameters to find the 
SUSY point. The technique has been pioneered in 7V=1 
SUSY with Wilson fermions by the DESY-Miinster group 
ficj |: here we discuss the extension to A/"=4 SYM. 

The supercurrent S^i is a linear combination of three 
dimension-7/2 operators. It is easy to find lattice op- 
erators which reproduce these continuum operators, at 
tree level and with contamination from higher dimension 
operators. The choices are not unique and at the non- 
perturbative level each lattice operator will mix with all 
continuum operators in the same symmetry channel. Dif- 
ferent choices of lattice operator will reproduce the con- 
tinuum operator with different normalization, mixings, 
and 0(a) suppressed higher dimension contamination. 
Hence we express the operators O^i in a continuum lan- 



guage, and leave the particulars of lattice transcription 
(which amounts to various "improvements" w.r.t. 0(a) 
discretization errors) for detailed studies. Our intention 
here is to lay out the methodology. 

In their analysis of the W=l SYM case, the DESY- 
Miinster group found two dimension-7/2 operators, the 
supercurrent and another fermionic current T^. These 
mix in the lattice-continuum matching and so one must 
write down two lattice operators with undetermined co- 
efficients in order to find something which corresponds 
purely to (plus 0(a) dimension-9/2 contamination). 
In the present 7V"=4 case there are 5 dimension-7/2 op- 
erators which we will name O^' 5 , and the renormalized 
7V=4 supercurrent will, in all generality, take the form: 

+{Z 4lv F^8 ij - Z 5 V2D^P L + ^P R ) A }^ A + 0{a) 
= Z n Ol i + 0(a), (9) 

where the terms on the righthand side are bare (lattice) 
operators. Note that: (1) at tree level the supercur- 
rent corresponds to Zi=Z2=Z^=l and Z4=Z^=Q; (2) 
the renormalization constants Z n are universal w.r.t. the 
index i due to the 5[/(4)^ symmetry preserved by the 
lattice. 

We can tune to the SUSY point by varying parameters 
to force correlation functions of this lattice-implemented 
d^S^^ to vanish up to O(a) corrections. Specifically, to 
tune two parameters we need to choose 6 operators 0™ i 
in the same symmetry channel as S^i (otherwise the cor- 
relation function vanishes automatically). The natural 
choice is O'-' plus one dimension-9/2 operator Q ^ i . 
One then measures the matrix of correlation functions 

M™"(i) = |d 3 f(O^(i,f)O^(0,0)) (10) 

whose t derivative is the correlation function between 
d ll O™ i and O n (0,i) at vanishing spatial momentum. 
Since the operators involved are dimension-7/2 we gener- 
ically expect the elements of M mn (t) to decay as i~ 7 . At 
the SUSY point and for the right choices of Z m , Z m M mn 
decays as ai~ 8 for all n. We can fix the undetermined ra- 
tios Z 2 ...5/Z1 by enforcing that this holds for n — 1 ... 4. 
Forcing that it hold for n=5, 6 gives two conditions which 
can be used to check whether we are at the SUSY point- 
tuning to the SUSY point is tuning for Z m M m5 ~ a t~ 8 
and Z m M me ^at~ 8 . Actually since one of the operators 
is dimension-9/2 we must force one linear combination 
Z m M mn c n to vanish as a 2 t~ 9 . 

We do not see an obstacle to using this procedure to 
tune more parameters, if it proves too difficult to tune 
some of the vanishing quartic couplings via the potential 
method. Therefore in principle the tuning to the SUSY 
point can be done by any mixture of the Ward identity 
method and the effective potential method. 
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We have seen that for 7V C =2,3 colors, there are four 
fine-tunings in the action. For N c >3 colors there are 
six. In addition, one must fix the four relative renor- 
malization constants in the supercurrent. All but one of 
the scalar potential counterterms can be fixed by match- 
ing the effective potential, as determined by the multi- 
canonical simulation, to the target theory scalar poten- 
tial Tr[0 m , <j) n ] [4> m i 4> n ] ■ The overall strength of this term 
cannot be determined from the effective potential, be- 
cause it will be expressed in terms of the bare operators 
in our approach. 

This leaves just six fine-tunings for all number of 
colors iV c : one fine-tuning of the bare kinetic coeffi- 
cient for the scalar, one overall scalar potential coef- 
ficient, and the four relative supercurrent coefficients. 
Thus a total of six Ward identities must be measured 
well enough to distinguish their simultaneous minimum 
w.r.t. Zfa Zi/Z 2 , ■ • • , Z4/Z5. 

Discussion 

By preserving the SU(A)r symmetry of the target the- 
ory, the number of counterterms that must be fine-tuned 
is greatly reduced. This can be done by implementing 
GW fermions, with the chirality of the Yukawa cou- 
plings implemented with auxiliary fermions, extending 
the method of [a]. Because counterterm fine-tuning can 
be isolated to the purely bosonic sector, it can all be 
done off-line, i.e., without the expense of fermion matrix 
inversions (the bottleneck for all dynamical fermion sim- 
ulations). This is a great advantage, because a very large 
number of points in the bare action parameter space will 
have to be scanned in order to find the N=A SYM point. 
Finally, we have explained how the overlap problem can 
be alleviated by taking a multicanonical approach, flat- 
tening the distributions that will be scanned over. 

The main limitation to this method is that, since A/"=4 
SYM is conformal, the continuum limit is not a weak cou- 
pling limit. Our proposal should work at weak coupling, 
where one knows that the infrared description will be 
in terms of the same degrees of freedom as one puts on 
the lattice. But there is no guarantee that one can find 
lattice parameters which correspond to strongly coupled 
continuum theories. 

The principal challenge is that the method re- 
quires GW fermions, which are numerically expensive- 
especially in a theory such as this one, with massless 
particles and the corresponding critical slowing down. It 
will be a challenge to generate enough configurations to 
measure quantities with sufficient accuracy to determine 
the SUSY point. Fermions are necessarily involved in 
the correlation functions of the supercurrent, so storage 
of propagators during the simulations will be essential to 
performing the fine-tuning w.r.t. Ward identities. The 



storage and computing resources that will be required 
will be substantial, but we believe that the exploratory 
studies that need to be done can be performed in the near 
term. For instance, one of the authors (JG) has access to 
the Computational Center for Nanotechnology Innova- 
tion (CCNI) at Rensselaer. Indeed, it is currently being 
used by Giedt and collaborators for M=l SYM simula- 
tions at a sustained actual compute rate of 1 Tflop/s, 
precisely the sort of resource that would be required to 
perform exploratory studies of this proposal. Obviously 
early stages of such work will be very much technical 
studies of the lattice theory. Continuum results will take 
much longer. Nevertheless, the beginnings of first princi- 
ples nonperturbative study of jV=4 SYM are not so far 
off, we believe, if the current proposal is pursued with 
some dedication and adequate resources. We hope to 
report on further progress in that direction in the near 
future. 
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